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The density-matrix renormalization group method is used to study the ground state of the two- 
chain zigzag-bond Hubbard model at quarter filling. We show that, with a proper choice of the signs 
of hopping integrals, the ring exchange mechanism yields ferromagnetic spin correlations between 
interchain neighboring sites, and produces the attractive interaction between electrons as well as 
the long-range pair correlations in the spin-triplet channel, thereby leading the system to triplet 
superconductivity. We argue that this novel mechanism may have possible relevance to observed 
superconductivity in Bechgaard salts. 
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More than two decades have passed since the discov- 
ery of superconductivity in so-called Bechgaard salts such 
as (TMTSF) 2 X with X=PF 6 , C10 4 , etc. Q, 0. The 
mechanism of superconductivity of this strongly corre- 
lated quasi-one-dimensional (ID) electron system [j| is, 
however, still an open issue. Recently, experimental evi- 
dences have been accumulating that the Cooper pairs of 
this system are in the spin-triplet state: they are from the 
measurements of the temperature dependence of Knight 
shift and spin-lattice relaxation ratio [J, 6] as well as 
that of the upper critical field 0, IM HI • 

A variety of theoretical approaches have so far been 
made on this unconventional superconductivity |lCl ITU 
IH fl3l fLU Hij . but none of them is based upon real-space 
pairing mechanism possibly relevant due to their strong 
electron correlations. In this paper, we want to propose 
a novel mechanism that may have possible relevance to 
triplet superconductivity of this system. We adopt a nu- 
merical approach in order to take fully into account the 
strong electron correlations. 

First, let us point out that the hopping integrals 



of (TMTSF)2X have the unique structure as shown in 
Fig. 1(a) [16|: they show an alternating sign change along 
the zigzag bonds connecting two chains, while the sign 
along the ID chain is always positive in the electron no- 
tation (i.e., for the bands of 3/4 filling of electrons). We 
hereafter use the hole notation for convenience; i.e., the 
system is in the quarter filling of holes. The signs of 
hopping integrals are then always negative along the ID 
chain and are changing signs along the zigzag bonds. 

Let us then notice that, in Hubbard models defined 
on such triangular-lattice related structures with proper 
signs of hopping integrals, the ring exchange mechanism 
of two spins on a triangle makes the system a ferromagnet 
with full spin polarization if sufficiently large Hubbard 
interaction U acts [l7l HI El Ej The structure of hop- 
ping integrals of (TMTSF^X |l6j. in fact, satisfies the 
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FIG. 1: Schematic representation of (a) the structure of hop- 
ping integrals for (TMTSF) 2 X and (b) the two-chain 
Hubbard model we use. The signs of hopping integrals are 
shown in the electron notation in (a) and hole notation in 
(b). The parameter values for (TMTSF) 2 C10 4 at low tem- 
peratures, e.g., are: tsi = 287, ts2 = 266, in = —34.0, 
t 12 = -64.1, and ti 3 = 46.2 meV 111. 



ferromagnetic sign rule txliH > (in the hole notation) 
for three hopping integrals of all the triangles. 

Our mechanism for triplet superconductivity relies on 
this ferromagnetic coupling: spin-triplet coupling for fer- 
romagnetism is 'relaxed' to short-range spin-triplet cor- 
relations in realistic strength of U, thereby leading the 
system to a metallic state with ferromagnetic spin corre- 
lations. We then ask what this state actually is. 

We calculate the ground state of the relevant two-chain 
Hubbard model by the density-matrix renormalization 
group (DMRG) method and show that indeed the 
system is metallic with ferromagnetic spin correlations. 
Furthermore, we show that the attractive interaction acts 
between holes, being caused by the gain in kinetic energy 
due to ring exchange of holes. We also show that the su- 
perconducting pair correlations in the spin-triplet chan- 
nel extend long-ranged in power-law length dependence, 



2 



while the pair correlations in the singlet channel as well 
as the spin-density-wave (SDW) correlations decay ex- 
ponentially, indicating that the system is in the state of 
spin-triplet superconductivity. We thus propose that the 
spin-triplet superconductivity is realized in the two-chain 
Hubbard model by the ring-exchange mechanism. Wc 
argue that the proposed mechanism may have possible 
relevance to triplet superconductivity in (TMTSF^X. 

The minimum model to include the effects of inter- 
chain coupling may be the two-chain Hubbard model (see 
Fig. 1(b)) defined by the Hamiltonian 



(1) 



<ij>a 



where c\ a (c Jcr ) is the creation (annihilation) operator of 
a hole with spin a (=|, j) at site i (j), = c\ a Ci a is the 
number operator, and (ij) denotes the nearest-neighbor 
pair. We restrict ourselves to the case at quarter fill- 
ing; i.e., n ~ ~^2 a { n ia) = 0.5 where (• • •} denotes the 
ground-state expectation value, tij is the hopping inte- 
gral between sites i and j: we include t\ along the zigzag 
chain and t 2 along the ID chains. U is the on-site Hub- 
bard repulsion. We also examine the intersite repulsive 
term J2<ij> Vij n i n j when necessary: we again include V\ 
along the zigzag chain and V 2 along the ID chains. We 
hereafter assume t 2 = — 1 as the unit of energy \t 2 \ = 1; 
observed small dimerization of hopping integrals along 
the ID chains is neglected because the system remains 
metallic when the interchain hopping integrals are of the 
zigzag type [i^. We use values t\ = ±0.25 and ±0.5 
with the sign alternation. The signs of t\ and t 2 can 
instead be taken all positive because the models where 
the product of the three hopping integrals of the tri- 
angles is positive are equivalent under canonical trans- 
formation. We use a value U = 10 (in some cases 20) 
for representing strong electron correlations in real sys- 
tems |23j . The DMRG method is used for calculating the 
ground-state energy and correlation functions for clusters 
of length L (containing 2L sites) with open boundary 
condition. We use clusters of up to L — 128 with keeping 
up to m ~ 4500 density-matrix eigenstates; the discarded 
weights are typically of the order 10~ 7 — 10~ 6 to obtain 
the ground-state energy in the accuracy of ^O.OOl^l- 
Lanczos exact-diagonalization method on small clusters 
is also used. 

Basic features of the ground state of our model are 
the following. First, the presence of ferromagnetic spin 
correlation is evident in Fig. 2, where we show the 
phase diagram on the (Vi,V 2 ) plane of an L=8 cluster 
with periodic boundary condition obtained by an exact- 
diagonalization method. We find that when U is large 
enough the system is spin polarized with the total spin 
S > (see Fig. 2(a)). When U is not large the total 
spin of the ground state becomes S — but the in- 
terchain nearest-neighbor spin correlation is still ferro- 
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FIG. 2: (a) Total-spin quantum numbers S of the ground 
state of the L — 8 cluster with U = 30 and = 1 plotted 
on the (Vi, V2) plane, (b) Contour plot of the spin correlation 
(Si • Sj) between the interchain nearest-neighbor sites for the 
L — 8 cluster with U — 10 and |ti| = 0.5 where Sj is the spin 
operator of a hole at site i. The ground state has S = 0. 



magnetic (see Fig. 2(b)); these results are confirmed to 
persist in larger size systems. The mechanism is appar- 
ent already in the two-hole ground state of a three-site 
Hubbard ring (which is spin triplet in a wide parameter 
region when t^t^ > 0) and persists in higher dimensions 
as well 0| . We also calculate the charge gap defined as 
A c = Hindoo A C (L) with A C (L) = E L (N + 2) + E L {N - 
2)—2El(N) where El(N) is the ground-state energy of a 
chain of length L with N electrons with equal number of 
t and i spins. We use the DMRG method for clusters of 
up to L — 64 and make a finite-size scaling analysis. The 
obtained results (not shown) are similar to the results 
obtained in Refs. [li], EH Ei| ; i.e., the charge gap opens 
due to charge ordering when either V\ or V 2 is large, but 
metallic state appears in a wide parameter region around 
Vi = 2V 2 (see Fig. 4 below for the case V\ = V2 = 0). 
Note that the region where the ferromagnetic correlation 
is strong agrees well with the metallic region where the 
ring exchange of holes works well. 

To see the behavior of the spin degrees of freedom fur- 
ther, we calculate the equal-time spin correlation func- 
tion S{i,j) = (SfS?) ~ (S?)(S?) where Sf is the z- 
component of the spin operator of a hole at site i. We also 
calculate the spin gap defined as A s = limi^ oc A s (L) 
with A S (L) = E L {N^ + l,i\T| - 1) - E^Nf,^) for 
2Vf = Ni = N/2. The results are shown in Fig. 3, where 
and hereafter we assume V\ — V 2 — unless otherwise 
stated. We first find that the values of A S (L) are rather 
small and extrapolated to A S (L) — > within our numeri- 
cal accuracy. Thus, the spin gap A s vanishes (or becomes 
quite small if it exists) in the thermodynamic limit. Wc 
also find that the spin correlation decays with nearly ex- 
ponential length dependence, which is associated with 
the oscillations of a period of 4 times lattice constant 
(consistent with the period of observed 2/cf-SDW state 
in (TMTSF) 2 X 0). We should note that the apparent 
contradiction between the vanishing spin gap and expo- 
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FIG. 3: Spin correlation function (a) at U = 10 and |ti| = 0.25 
and (b) at U = 10 and |ti| = 0.5 for the L = 128 cluster. 
Horizontal axis \i — j\ counts the sites along the zigzag chain, 
(c) Spin gap A S (L) at U — 10 plotted as a function of 1/L. 
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FIG. 4: Binding energy A b (L) and charge gap A c (L)/4 plot- 
ted as a function of 1/L. (a) U — 20 and ti = 0.5, and 
(b) U = 10 and = 0.5. Solid curve connecting A C (L) is 
the fitting by a polynomial of 1/L. Solid and dotted lines 
connecting Aj~ (L) are a guide to the eyes. 



nential decay of the spin correlation (as well as its rapid 
decrease at \i — j\ < 30) may be reconciled by the spin- 
triplet pairing present in our model (see below) [2t| • We 
point out that the Tomonaga-Luttinger-liquid descrip- 
tion [28| of our results obtained in the situation where 
there are four Fermi poi nts in the strong coupling regime 
is yet an open issue [l9j . 

We now calculate the binding energy of holes defined 
as A fc = Hindoo A±(L) with Af(L) = E L (N ± 2) + 
E L (N) - 2E L (N ± 1). Because E L (N) depends only on 
n = N/ (2L) for a given system of length L, we have the 
scaling property for large L as 



A± 



(L) 



1 d 2 (E L /L) 
4L dn 2 



L 3 ' 



(2) 
(3) 



where the proportionality coefficient of 1 /L is related to 
the inverse of compressibility n as 



= n lim 



d 2 (E L /L) 
dn 2 



(4) 



The binding energy A& in the thermodynamic limit is 
obtained in Eq. (3) as the 0(1) correction in Af(L) [2^| . 
The O(l) correction in A C (L), i.e., A c , vanishes when 
the system is metallic. 

The calculated results are shown in Fig. 4. We first 
note that the values of k evaluated from the calculations 
of A C (L) and Af(L) agree well with each other as the 
gradients of the two curves (if the factor 4 is taken into 
account) for small 1/L regions agree well. Also noted is 
that the gradient is positive, or k > 0, indicating that the 



system is thermodynamically stable against phase sep- 
aration. We then find that the extrapolated value A;, 
is negative; i.e., the attractive interaction acts between 
holes in the thermodynamic limit. The energy gain re- 
sponsible for the negative value of A& may come from the 
motion of two holes running around the triangle, avoiding 
the on-site repulsion U and exchanging spins for triplet 
coupling, i.e., from the ring exchange of holes. The pair- 
ing mechanism is thus kinetic in origin. 

We also calculate the pair correlation function defined 
as D(l) = (At + jAi) with A, = c^c i+r i - Ciic i+r ^ for sin- 
glet pairs and Aj = c^Ci+ r | + c^Ci+ r | for triplet pairs 
where i + r denotes the neighboring sites of i. The re- 
sults at L — 128 are shown in Fig. 5. We find that D(l) 
shows the power-law length dependence for the interchain 
triplet pairing but decays exponentially for the singlet 
pairing as well as for the triplet pairing on the single ID 
chain. These are the case also at L = 64, indicating that 
the size of the clusters used is sufficiently large. Quan- 
titatively, the pair correlation function (see Fig. 5(b)) at 
long distances decays as ~ 1/rj- 8 , the decay of which 
is slower by comparison than the decay of the charge 
correlation function ~ l/r 2 j°; the estimation is made by 
fitting the data at 1 <C I <C L in order to avoid the ef- 
fects of the edges of the clusters. Then, combined with 
the attractive interactions of two holes shown above, our 
results indicate that the system should be in the state of 
spin-triplet superconductivity where the pairing of holes 
occurs between the interchain nearest-neighbor sites. 

Finally, let us consider possible relevance of our results 
to Bechgaard salts. The model is a coupled sequence of 
the ID chains. We have cut out two neighboring chains 
as a minimum model to seek for consequences of the in- 
terchain coupling. Since the ring-exchange mechanism 
works also for models of more than two chains 181 where 
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FIG. 5: (Color online) Pair correlation functions D(l) calcu- 
lated at L = 128 and U = 10 with (a) \ti\ = 0.25 and (b) 
\t\\ = 0.5. Horizontal axis I counts the sites along the ID 
chain. In the inset, the two circles (squares) denotes the sites 
i and i + I (i + r and i + l + r), where 1=1 corresponds to the 
position for two squares to sit on the same site. 



the hopping integrals of all the triangles satisfy the ferro- 
magnetic sign rule, we point out that the triplet pairing 
state obtained in the two-chain model may persist also in 
quasi-lD (or 2D) systems; future studies will be interest- 
ing. Also pointed out is that there is intriguing competi- 
tion between the SDW and spin-triplet superconducting 
states in the experimental pressure-temperature phase di- 
agram [3(J ■ Because small intersite Coulombic repulsions 
exist in the materials |23j|. we examine their effects; we 
find that the inclusion of a realistic value V% ~ 1 
does not change our results. However, if we include V\ 
between the ID chains as well, the triplet pair correlation 
becomes less long-ranged and thus the SDW correlation 
can be comparable with it. A true long-range order may 
be selected among these competing correlations when the 
two (or three) dimensionality is taken into account, al- 
though to predict which order is realized is beyond the 
scope of the present work. We hope that future quanti- 
tative analyses will help clarifying the issue. 

In summary, we have calculated the binding energy 
of holes and pair correlation functions in the two-chain 
zigzag-bond Hubbard model and have shown that the 
model can be superconducting in spin-triplet channel 



when we make an appropriate choice of the signs of hop- 
ping integrals for the ferromagnetic ring-exchange mech- 
anism to work. We have argued that the mechanism 
proposed may have possible relevance to the triplet su- 
perconductivity in (TMTSF) 2 X. 
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